function [bo] = create_hematite_rhombohedral()
%CREATE_HEMATITE_RHOMBOHEDRAL Summary of this function goes here
%   Detailed explanation goes here
z = 0.10534;
x = 0.3056;
basisVectors = [z z z;
  1/2-z 1/2-z 1/2-z;
  1/2+z 1/2+z 1/2+z;
  1-z   1-z 1-z;
  x 1-x 0;
  1-x 0 x;
  0 x 1-x;
  1/2-x 1/2+x 1/2;
  1/2+x 1/2 1/2-x;
  1/2 1/2-x 1/2+x];
atoms = {'Fe';
  'Fe';
  'Fe';
  'Fe';
  'O';
  'O';
  'O';
  'O';
  'O';
  'O';};
bo = Basis(basisVectors, atoms);

end

